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Abstract 

Polarization-resolved optical spectra of magnesia-supported gold clusters AuTv/MgO 
(N=l, 2,4,8), bound at a surface color center F s of the MgO(lOO) face, are calculated from the 
time-dependent density functional theory. The optical lines for N=l,2 are dominated by transi- 
tions that involve strong hybridization between gold and F s states whereas for iV=4,8 intracluster 
transitions dominate. The theoretical optical spectra are sensitive to cluster structure and adsor- 
bants (here CO and O2 molecules on Aus/i^QMgO) which suggests polarization-resolved optical 
spectroscopy as a powerful tool to investigate structures and functions of chemically active, sup- 
ported clusters. 

PACS numbers: 36.40.Vz,68.47.Jn,36.40.Mr,73.22.Lp 
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I. INTRODUCTION 



Investigations of the physical and chemical properties of metal clusters are currently 
largely motivated by the question how their various remarkably size- dependent properties 
could be best utilized while the clusters are interacting with the environment, e.g., bound on 
or implanted in a support, or stabilized and surface-passivated by molecules^. Understanding 
factors that dictate the stability, structure and function has relevance regarding atomic-scale 
design of components that could be of potential use in future nanotechnologies. To this end, 
spectroscopic tools and density functional theory (DFT) calculations can provide valuable 
insights. 

Gold clusters and nanoparticles have attracted much attention recently due to their re- 
markable chemical and optical properties (for recent reviews, see e.g.-^). Gold clusters with 
a size of just a few atoms, supported by thin MgO films, have turned out to be paradigms to 
provide insight into the importance of quantum finite size effects in nanoscale chemistry - it 
has been demonstrated that by changing the size or elemental composition of the supported 
cluster even atom by atom one can quite dramatically affect the catalytic properties of the 
cluster— Earlier DFT calculations have shown that the finite-size effects in the electronic 
spectrum of the nanocatalyst play a key role in CO oxidation reaction through binding and 
activating the oxygen molecule via a charge-transfer mechanism. Theoretical predictions^ 
of charging of the catalytically active Au 8 cluster, bound at the surface color center F s of 
the MgO film, were confirmed very recently by FTIR spectroscopy^. 

Spectroscopic data on structure of supported metal clusters and ultrafine nanoparticles 
are still relatively scarce, with some notable exceptions from STM and STS spectroscopy^. 
Optical spectroscopy is currently a largely unexplored area, although it has long been used to 
study metal clusters in gas phase*- Here we present a systematic theoretical investigation of 
the optical properties of MgO-supported gold clusters, Aujv/MgO with N < 8, by analyzing 
the optical spectra calculated from the linear response time-depenent DFT (TDDFT). We 
discuss (i) the effects from the support, particularly from the F s states, dominating the spec- 
tra of N = 1,2, (ii) the sensitivity of the spectra to the atomic structure of the cluster with 
examples for N = 4, 8, and finally (iii) the sensitivity of the spectra to adsorbants, by taking 
the example of CO and O2 adsorption on the smallest CO-oxidizing cluster Aus/Mg Q 4 ' 5 ' 6 . 
It is shown that for all these cases, polarization-resolved spectra contain rich information 
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on the nature of the dominant optical transitions and their connection to the structure and 
shape of the adsorbed cluster. 

II. GROUND-STATE DFT CALCULATIONS 

The atomic and electronic structure of the supported Au/v/MgO system (comprising the 
vicinity of the color center of the MgO(lOO) surface, the adsorbed Au^r cluster, and O2 and 
CO molecules) were calculated within the DFT in combination with Born-Oppenheimer 
(BO) Molecular Dynamics (MD)^ including self-consistent gradient corrections via the so- 
called PBE-GGA functional 11 . Au(5rf 10 6s 1 ), Mg(3s 2 ), C(2s 2 2p 2 ), and 0(2s 2 2p 4 ) electrons 
were included in the valence, and the interaction to the ion cores was described by scalar- 
relativistic non-local norm conserving pseudopotentials devised by Troullier and Martin a 12 ' 13 . 

The MgO surface is modelled by a two-layer ab initio cluster Mg m O m or F s @Mg m O m _i, 
embedded in an extended point-charge lattice to include effects of the long-rage Madelung 
potential. For the embedding lattice, around 2100 alternating charges of +2 and -2 rep- 
resenting Mg and O ions, respectively, were used. In addition, those positive point charges 
that would be nearest neighbors to the peripheral O atoms of the central ab initio Mg m O n 
cluster have been replaced by "empty" Mg pseudopotentials (MgPP) in order to prevent un- 
physical polarization effects to O ions (see visualization of Au/Mg 13 0i3 in Fig. la) 14 . The 
lattice parameter of the embedding part is fixed to the experimental lattice constant (4.21 
A) of bulk MgO. The cluster, molecules and nearest-neighbor Mg ions to the F s are treated 
dynamically in structural optimizations that included both steep est- descent and quenched 
molecular dynamics runs^* 1 ^ 

III. LINEAR-RESPONSE TDDFT CALCULATIONS 

We follow the formulation of the linear response time dependent DFT given by Casida/ 6 
as implemented in ref.— . Briefly, in order to get the weights Fj and energies fiwi of optical 
transitions {/}, one solves an eigenvalue problem QFi = ujFi where the Q matrix elements 
are given by 

&ij,kl = Sik&jlt'ij + ly/nijEijUklEklKij^kl- (1) 
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Here = Ej — £j and ntj = n» — are the difference of the KS particle-hole eigenvalues 
and occupation numbers, respectively. K^^i is a coupling matrix that describes the linear 
response of the electron density p to the single-particle - single-hole excitations in the basis 
spanned by the ground state KS orbitals The transition matrix element in polarisation 
direction i v is 



and the corresponding polarization-dependent and polarization-averaged oscillator strengths 
are (//)„ = 2\(M I ) U \ 2 and // = (1/3) Y^v=i(fi)vi respectively. The convergence of the spectra 
(oscillator strengths and energies of the major peaks) depends basically on three factors: (i) 
extent of the KS particle-hole basis, (ii) size of the plane-wave grid, and (iii) size of the 
ab initio Mg m O n cluster. Factors (i) and (ii) determine the quality of discretization of 
the continuum virtual states while factor (iii) dictates the overall behavior in the optical 
response of the MgO support - specifically the width of the band gap. We have carefully 
tested the convergence behavior of the calculated spectra in case of the smallest supported 
clusters (N—1,2) by using Mg m O m _i clusters extending up to m = 49. The converged 
theoretical optical band gap of the support is about 4 eV, which is in a fair agreement with 
the apparent gap of 6 eV measured for a few- layer thick MgO film 18 , taking into account 
the general tendency of DFT to underestimate band gaps of wide-gap insulators. However, 
the important transitions induced by the supported cluster fall into the band gap. These 
transitions can be calculated within 0.1 eV accuracy even with significantly smaller cluster 
sizes (m = 13, 25). 

IV. RESULTS AND DISCUSSION 

Application of the TDLDA method to an isolated gold atom leads to the optical spectrum 
shown in the inset to Fig. lb (for comparison, we show also the optical spectrum based on 
KS-only transitions, i.e., by setting K^^i = in Eq. (1)). This spectrum was calculated by 
considering the (5d l0 6s 1 6p°7 s°7p°) KS states in the particle-hole basis. The KS spectrum 
consists of two strong transitions corresponding to single-particle excitations of 6s — > 6p 
(5 eV) and 5d — > 6p (6.1 eV). The lowest TDLDA optical line is found at 5.23 eV and it 
is an incoherent (destructive) superposition of the KS transitions 6s — > 6p and 5c? — > 6p 
(the destructive interference is manifested by the reduced oscillator strength of the 5.23 
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FIG. 1: (color online) The supported Au atom: a) The ab initio cluster Au/Mgi30i3 with 
the peripheral "empty" MgPP sites, b) the polarization averaged spectrum (the inset shows the 
spectrum of a gas phase Au atom), c) the polarization resolved spectrum, d) and e) are polarization 
averaged and resolved spectra for the atom on F s . All spectra here and in subsequent figures are 
plotted by using a Lorenzian folding of 0.08 eV of individual optical lines. Note, that the spectra 
in d,e have been calculated for Au/F s @Mg2s024. 

eV TDLDA line compared to the neighboring KS lines). Our result of 5.23 eV is in a 
fair agreement with the lowest experimental spin-orbit (SO) split transitions at 4.63 eV 
( 2 Si/2 — > 2 Pi/2, final single-particle state assigned as (5d w 6p 1 )) and 5.11 eV ( 2 Si/ 2 — > 
2 Pz/2 (5<i 10 6p 1 ))— taking into account the fact that the SO effect is not included in our 
calculations. Similarly, the next strong TDLDA transition at 7.3 eV corresponds to the 
experimental SO doublet 7.44 eV and 7.53 eV, assigned to transitions to the 2 -Pi/2, 2 P^/2 
states with a (5<i 10 7p 1 ) single-particle configuration 19 . We also note here, that another recent 
TDDFT calculation 2 ^ finds the first optical transition of an isolated Au atom at 5.23 eV, in 
full agreement with our result. 

For the combined Au/MgO system, we find the lowest strong transitions at 3.38 eV and 
3.61 eV, i.e, in the near-UV region (see Fig. lb). The nature of these lines is revealed 
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FIG. 2: (color online) a) Two configurations for Au2 bound to F s , b) polarization averaged and 
c,d) resolved spectra. d(Au-Au)=2.61 A and 2.62 A for "up" and "tilted" geometries, respectively, 
compared to the gas-phase value of 2.55 A. 

by the the polarization-resolved spectra (Fig. lc) which shows basically three peaks that 
have degenerate (fi) v values, two of which (x,y polarization) at 3.38 eV and the z-polarized 
transition (i.e., excitable for z-polarized light) at 3.61 eV. The identical oscillator strengths 
in the x,y,z directions and the degeneracy pattern strongly suggest that these transitions 
mainly originate from the Au atom, perturbed and shifted by the support. (Note that 
neither the ideal MgO surface nor the gas phase atom have transitions anywhere near 3.5 
eV). One can also locate a set of weak transitions in the region 1.7 eV - 1.9 eV that consists 
of three weak lines, of which two are z-polarized and one xy-polarized. In general, these 
lines are due to incoherent mixing of several KS transitions, but we note that they have 
a component of the atomic hd — > 6s character, which is dipole-forbidden in the gas phase 
atom but becomes partially "visible" via mixing with the substrate states. 

The situation is dramatically changed when the F s defect is introduced (Fig. ld,e). 
Several strong lines are now seen in the range 1.8 eV - 3.5 eV. Fig. le reveals strikingly 
that basically all of the transitions are polarized only in the z-direction, i.e., perpendicular 
to the surface. A stronger coupling to the F s states also enhances the oscillator strength of 
the bands around 2 eV which have the contribution from the gas-phase forbidden 5c? — > 6s 
transition. 

For Au 2 , bound at F s , we have found two almost isoenergetic geometries (Fig. 2a); the 
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FIG. 3: (color online) The configuration a) and polarization resolved spectrum b) of AU4 bound 



"tilted" geometry with the Au-Au axis 33° from the surface normal is 0.1 eV more stable 
than the dimer standing "up". This slight difference in the adsorption geometry shows up 
drastically in the TDLDA spectrum (Figs. 2b-d): the "up" geometry has transitions only in 
the z-direction whereas the spectrum of the "tilted" cluster has a clear y-component at 3.1 
eV (y is the direction to which the dimer is tilted). Based on the clear difference between the 
TDLDA spectra in Figs. 2b-d we predict that in this simplest non-trivial case of adsorption 
of a gold cluster on the defect site of MgO surface, experimental optical spectra could yield 
a definitive answer on the binding geometry. 

AU4, whose optimum binding geometry has previously been found to be that of a tilted 
rhombus-like structure^ (Fig. 3a), has a very strong transition close to 2.5 eV (Fig. 3b). 
This transition is solely in the direction of the long axis of the cluster. The tetramer is thus 
the smallest cluster which "screens out" the transitions involving major (z-) components from 
the color center, which had been found to be dominant for the monomer and dimer. 

We now turn to discussion of the optical properties of the "smallest MgO-supported 
nanocatalyst" for CO oxidation, namely Aus-^^ Two isomeric structures were discussed 
in Ref. 5 , where the open structure A is more stable than the compact structure B by 0.3 
eV (see Fig. 4a, note that similar geometric motifs yield locally stable isomers also for 
neutral gas-phase gold octamers where open structures are favored due to relativistic bonding 
effect o 22123 ). Fig. 4b shows that these two isomers yield very different, characteristic spectra: 
while the more stable open structure A has a series of intense transitions covering most of 
the optical range and yielding characteristic line shapes, the spectrum of structure B has less 
oscillator strength below 2.7 eV. Fig. 4c gives more insight into the strong absorption bands 
at 2.1 eV, 2.4 eV, and 2.65 eV of structure A: they are mostly transitions involving the long 
axis (~y-direction) of the cluster. Bands after 3 eV start to have major contributions also 
from transitions in the other two polarization directions. Again, Figs. 4b-c suggests that in 
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FIG. 4: (color online) Aus on F s : a) Two configurations and b) their polarization averaged spectra, 
c) The polarisation resolved spectra of structure A. 
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FIG. 5: (color online) Aus and absorbed molecules on F s : a) The configurations of the adsorbed 
molecules and the spectra of b) structure A with and without O2 and c) structure B with and 
without CO/0 2 . 

principle it should be possible to distinguish different isomers of adsorbed gold octamers via 
high-resolution optical data. 

Finally, we show examples about modification of the optical spectra by adsorbant 
molecules in Fig. 5, here CO and O2. Fig. 5b shows that adsorption of the O2 molecule 
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TABLE I: Local charges^ (in units of e) on the constituents of the structures shown in Figs. 4a) 
and 5a). Support includes the MgO surface with the color center. 

to the structure A of the gold octamer strongly suppresses the absorption bands at 2.1 eV, 
2.4 eV, and 2.65 eV, furthermore, a new clear absorption band appears now at 1.5 - 1.7 
eV. Analysis of the polarization-resolved spectra (not displayed) shows disappearance of the 
bands connected to excitations along the long axis of structure A. On the other hand, ad- 
sorption of CO or 2 to the structure B of Aug does not significantly modify the absorption 
spectrum below 3 eV (Fig. 5c). This finding can be correlated to charge transfer horn Aug 
to 2 (see Table 1): local charge on the MgO-supported gold octamer (structure A) changes 
from — 0.37e to +0.73e upon oxygen adsorption, and 2 gains a charge of — 1.26e that goes 
to the originally empty 2tt* antibonding orbitals, i.e., the adsorbed dioxygen is activated to 
a peroxo-like state.— This charge-transfer-induced activation was shown to be the key to 
the catalytic acitivity.^^ In the TDDFT framework, the significant change of the oxidation 
state of the gold octamer A upon dioxygen adsorption (Lie in total) obviously depletes 
those single-particle KS states from the TDDFT basis, that have a major contribution to 
the strong excitations of the bare cluster in the energy range shown. 2 and CO adsorption 
on structure B are accompanied by a smaller change of the oxidation state of gold (0.8e and 
0.58e in total, respectively). 

We have shown that polarization-resolved optical spectra of magnesia-supported gold 
clusters contain rich information about the binding mode and structure of the cluster and 
possible adsorbant molecules. Specifically, optical spectra contain complementary infor- 
mation about the adsorption of CO and 2 molecules on the catalytically active center 
Aug/F s @MgO and the related charge transfer. Emerging applications of surface-sensitive 
optical spectroscopic tools 20 to this and other chemically active supported metal clusters are 
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likely to help significantly in gaining detailed understanding of structures and functions of 
reaction centers of nano catalysts. 
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